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Abstract. Envariance - entanglement assisted invariance - is a recently discovered 
symmetry of composite quantum systems. We show that thermodynamic equilibrium 
states are fully characterized by their envariance. In particular, the microcanonical 
equilibrium of a system S with Hamiltonian H$ is a fully energetically degenerate 
quantum state envariant under every unitary transformation. The representation of 
the canonical equilibrium then follows from simply counting degenerate energy states. 
Our conceptually novel approach is free of mathematically ambiguous notions such as 
ensemble, randomness, etc., and, while it does not even rely on probability, it helps to 
understand its role in the quantum world. 


1. Introduction 

Statistical physics was developed in the XIX century. The fundamental physical theory 
was then Newtonian mechanics. The key task of statistical physics was to bridge the 
chasm between microstates (points in phase space) and thermodynamic macrostates 
(given by temperature, entropy, pressure, etc.) Strictly speaking, that chasm is - 
within the context of Newtonian mechanics - unbridgeable, as classical microstates 
have vanishing entropy. Thus, a “half-way house” populated with fictitious but useful 
concepts such as ensembles was erected - half way between micro and macro - and 
served as a pillar supporting the bridge. Even before ensembles were officially introduced 
by Gibbs [1], the concept was de facto used by, e.g., Maxwell [2] and Boltzmann [3]. 
Doubts about this “half-way house” strategy nevertheless remained, as controversies 
surrounding the H-theorem demonstrate. 

Our point is that, while in XIX century physics ensembles were necessary because 
thermodynamics anticipated the role played by information in quantum physics, the 
state of a single quantum system can be mixed. Therefore, the contradiction between the 
pure classical microstate and an (impure) macrostate does not arise in a Universe that 
is quantum to the core. Yet, the development of quantum statistical physics consisted 
to a large extent of re-deploying strategies developed to deal with the fundamental 
contradiction between Newtonian physics and thermodynamics. 
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We claim that such “crutches” that were devised (and helpful) in the XIX century 
became unnecessary with the advent of quantum physics in the XX century. Thus, in 
the XXI century we can simply dispose of the ensembles invented to justify the use of 
probabilities representing the ignorance of the observers and to compute the entropy of 
the macrostate. 

In the present paper we propose an alternative approach to the foundations of 
statistical mechanics that is free of the conceptual caveats of classical theory, and 
relies purely on quantum mechanical notions such as entanglement. Our approach is 
based on envariance (or entanglement assisted invariance) - a symmetry-based view 
of probabilities that has been recently developed to derive Born’s rule from the non- 
controversial quantum postulates [4-8] - and that is, therefore, exceptionally well-suited 
to analyze probabilistic notions in quantum theories. 

1.1. Entanglement assisted invariance 

Consider a quantum system, S , which is maximally entangled with an environment, £, 
and let | ifss) denote the composite state in 5®^. Then \ifss) is called envariant under 
a unitary map Us = us <S> I s iff there exists another unitary Us — Is ® us such that, 

Us | f’ss) — {us <S> If) | ifss) — I Vss ); (1) 

Us | rjss) = (Is ® us) \vss) = | Dss) ■ (2) 

Thus, Us that does not act on S “does the job” of the inverse map of Us on S - assisted 
by the environment £. 

The principle is most easily illustrated with a simple example. Suppose S and £ 
are each given by two-level systems, where {IDs, IDs} are the eigenstates of S and 
{lt)f ) IDs! span £. Now, further assume \ifse) oc lT)s ® IDs + IDs ® IDs an( 4 Us is a 
swap in S - it “flips” its spin. Then, we have 


IDs®IDf +IDs®IDs IDs ® IDs + IDs ® IDs 


( 3 ) 


The action of Us on | if) ss can be restored by a swap, I/s, on £, 


IDs^lDs+IDs^lDs ^ IDs ® IDs + IDs® IDs■ 


( 4 ) 


Thus, the swap Us on £ restores the pre-swap \if) S£ without “touching” S, i.e., the 
global state is restored by solely acting on £. Consequently, local probabilities of the 
two swapped spin states are both exchanged and unchanged. Hence, they have to be 
equal. This provides the fundamental connection of quantum states and probabilities [6], 
and leads to Born’s rule. 
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Envariance of pure states is a purely quantum symmetry. Classical pure states of 
composite systems are given by Cartesian rather then tensor products. Therefore, it 
can be shown that such an environment assisted inverse of a map acting on a classical 
system cannot exist [6]. This reflects the fact that some properties of quantum states 
can be considered “relative” with respect to S, whereas all properties of classical states 
are always “absolute”. 

Recent experiments [9, 10] have shown that envariance is not only a theoretical 
concept, but a physical reality. Thus, in the present work we will use envariance to 
advance a conceptually novel approach to the foundations of statistical mechanics. 

2. Microcanonical state from envariance 

We begin by considering the microcanonical equilibrium. Generally, thermodynamic 
equilibrium states are characterized by extrema of physical properties, such as maximal 
phase space volume, maximal thermodynamic entropy, or maximal randomness [11]. 
We will define the microcanonical equilibrium as the quantum state that is “maximally 
envariant”, i.e., envariant under all unitary operations on S. To this end, we write the 
composite state | ipse) hi Schmidt decomposition, 

I ipse) = ak l Sfc ) ® ’ ( 5 ) 

k 

where by definition {|>Sfc)} and (|£fc)} are orthocomplete in S and £, respectively. The 
task is now to identify the “special” state that is maximally envariant. 

It has been shown [6] that \ifss) is envariant under all unitary operations if and 
only if the Schmidt decomposition is even, i.e., all coefficients have the same absolute 
value, |cifc| = |a/| for all l and k. We then can write, 

I'f’se) oc ^2 ex P {ifa) |sfc) ® kfc) , (6) 

k 

where </>*. are phases. Recall that in classical statistical mechanics equilibrium ensembles 
are identified as the states with the largest corresponding volume in phase space [11], 
In the present context this “identification” readily translates into an equilibrium state 
that is envariant under the maximal number of, i.e., all unitary operations. 

It is interesting to note that Eq. (6) can also be used to derive Born’s rule [4,5]. In 
particular, it can be shown that an outside observer would measure each of the states 
| Sk) with equal probability. However, we emphasize once again that the present analysis 
is independent of any notion of probability, and all steps rely purely on symmetries of 
entanglement. 
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To conclude the derivation we note that the microcanonical state is commonly 
identified as the state that is also fully energetically degenerate! [12]. To this end, 
denote the Hamiltonian of the composite system by 

Hse — Hs ® If + Is ® He ■ (7) 

Then, the internal energy of S is given by the quantum mechanical average 

E s = {il>se\ (H s <8> h) 1 1>se) = < s *l H s k> ! z > (8) 

k 

where Z is the energy-dependent dimension of the Hilbert space of S , which is commonly 
also called the microcanonical partition function [12]. Since | f’ss) (6) is envariant 
under all unitary maps we can assume without loss of generality! that {sk} f =1 is a 
representation of the energy eigenbasis corresponding to H s , and we have (s*, | H s |s fc ) = 
ek with Es = ek = for all k, k' G {1,..., Z}. 

In conclusion, we have identified the fully quantum mechanical representation of 
the microcanonical state by two conditions§: the state representing the microcanonical 
equilibrium of a system S with Hamiltonian Hs is the state that is (i) envariant under 
all unitary operations on S and (ii) fully energetically degenerate with respect to H s . 

3. Reformulation of the fundamental statement 

Before we continue to rebuild the foundations of statistical mechanics using envariance, 
let us briefly summarize and highlight what we have achieved so far. All previous 
treatments of the microcanonical state relied on notions such as probability, ergodicity, 
ensemble, randomness, indifference, etc. However, in the context of statistical 
physics none of these expressions are fully well-defined. Indeed, in the early days of 
statistical physics seminal researchers such as Maxwell and Boltzmann struggled with 
the conceptual difficulties [11]. Modern interpretation and understanding of statistical 
mechanics, however, was invented by Gibbs, who simply ignored such foundational issues 
and made full use of the concept of probability. 

In contrast, in our approach we only need a quantum symmetry induced by 
entanglement - envariance - instead of relying on mathematically ambiguous concepts. 
Thus, we can reformulate the fundamental statement of statistical mechanics in quantum 
physics: 

The microcanonical equilibrium of a system S with Hamiltonian Hs is a fully 
energetically degenerate quantum state envariant under all unitaries. 

The remainder of this analysis will further illustrate this novel conceptual approach to 
the foundations of statistical mechanics by also treating the canonical equilibrium. 

! We would like to emphasize that this is a standard assumption in textbook derivations [12]. It can 
be easily relaxed (also employing standard tools) by considering an energy interval [13]. For related 
discussions in a context similar to the present work see also Refs. [14,15]. 

§ Note that in our framework the microcanonical equilibrium is not represented by a unique state, but 
rather by an equivalence class of all maximally envariant states with the same energy. 
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4. Canonical state from quantum envariance 

Let us now imagine that we can separate the total system S into a smaller subsystem 
of interest © and its complement, which we call heat bath 23. The Hamiltonian of S 
can then be written as 


H$ — Hq © I® + Ig <S> H<% + he,® , (9) 

where denotes an interaction term. Physically this term is necessary to facilitate 
exchange of energy between the © and the heat bath 23. In the following, however, 
we will assume that h©,® is sufficiently small so that we can neglect its contribution 
to the total energy, E s = E & + E<g, and its effect on the composite equilibrium state 
| ifss)- These assumptions are in complete analogy to the ones of classical statistical 
mechanics [12,13]. They will, however, be relaxed in a final part of the analysis. 

Under these assumptions every composite energy eigenstate | Sk) can be written as 
a product, 

|sfc) = |fifc) <8 |b fc ) , (10) 

where the states |Sfc) and |bfc) are energy eigenstates in © and 23, respectively. At this 
point envariance is crucial in our treatment: All orthonormal bases are equivalent under 
envariance§. Therefore, we can choose |sfc) as energy eigenstates of Hg. 

For the canonical formalism we are now interested in the number of states accessible 
to the total system S under the condition that the total internal energy Es (8) is given 
and constant. When the subsystem of interest, ©, happens to be in a particular energy 
eigenstate |sfc) then the internal energy of subsystem is given by the corresponding 
energy eigenvalue e*,. Therefore, for the total energy Es to be constant, the energy of 
the heat bath, E<g, has to obey, 

E<& ( e fc) = E s — C/c • (11) 


This condition can only be met if the energy spectrum of the heat reservoir is at least 
as dense as the one of the subsystem. This observation will be illustrated shortly in 
Eqs. (14) and (16). 

The number of states, 21 (c*,), accessible to S is then given by the fraction 


21 (cfc) 


21 ® (Es — tk) 
21 s(E s ) 


( 12 ) 


where 21 s(Es) is the total number of states in S consistent with Eq. (8), and 
21*8 (Eg — Cfc) is the number of states available to the heat bath, 23, determined by 
condition (11). In other words, we are asking for nothing else but the degeneracy in 23 
corresponding to a particular energy state of the system of interest |5fc). 


4-1. Example: Composition of multiple qubits 

The idea is most easily illustrated with a simple example, before we will derive 
the general formula in the following paragraph. Imagine a system of interest, ©, 
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that interacts with N non-interacting qubits with energy eigenstates |0) and |1) and 
corresponding eigenenergies ef and ef. Note once again that the composite states | s k ) 
can be always chosen to be energy eigenstates, since the even composite state \ifss) (6) 
is envariant under all unitary operations on <S§. 

We further assume the qubits to be non-interacting. Therefore, all energy 
eigenstates can be written in the form 

I s k ) = | s k ) <g> j • (13) 

N—qubits 

Here 8 l k G {0,1} for all i G 1,..., N describing the states of the bath qubits. Let us 
denote the number of qubits of 53 in |0) by n. Then the total internal energy E$ becomes 
a simple function of n and is given by, 

E s = t k + nef + (N-n)ef. (14) 


Now it is easy to see that the total number of states corresponding to a particular value 
of Es, i.e., the degeneracy in 53 corresponding to e*,, (12) is given by, 


51 (e fc ) = 


N\ 


n\ (N — n)\ 


(15) 


Equation (15) describes nothing else but the number of possibilities to distribute nef 
and (N — n) ef over N qubits. 

It is worth emphasizing that in the arguments leading to Eq. (15) we explicitly used 
that the | s k ) are energy eigenstates in S and the subsystem G and heat reservoir 53 are 
non-interacting (10). The first condition is not an assumption, since the composite \ifss) 
is envariant under all unitary maps on S, and the second condition is in full agreement 
with conventional assumptions of thermodynamics. 


4-2. Boltzmann’s formula for the canonical state 

The example treated in the preceding section can be easily generalized. We again 
assume that the heat reservoir 53 consists of N non-interacting subsystems with identical 
eigenvalue spectra {e®}”L 1 . This model includes the most relevant cases commonly 
analyzed in quantum thermodynamics, such as quantum Brownian motion [16], where 
53 consists of N non-interacting harmonic oscillators. 

In this case the internal energy (8) takes the form 

E s = tk + ni ef + n 2 ef H-b n m e® , (16) 

with Xq=i n j = N- Therefore, the degeneracy (12) becomes 

51 (cfc) = 


ni!n 2 ! • ■ ■ n 


(17) 
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This expression is readily recognized as a quantum envariant formulation of Boltzmann’s 
counting formula for the number of classical microstates [11], which quantifies the volume 
of phase space occupied by the thermodynamic state. However, instead of having to 
equip phase space with an (artificial) equispaced grid, we simply count degenerate states. 

We are now ready to derive the Boltzmann-Gibbs formula. To this end consider 
that in the limit of very large, N 1, 91 (e*,) (17) can be approximated with Stirling’s 
formula. We have 

m 

In (91 (c fc )) = N In (N) — nj In (nf) . (18) 

3 = 1 

As pointed out earlier, thermodynamic equilibrium states are characterized by a 
maximum of symmetry or maximal number of “involved energy states”, which 
corresponds classically to a maximal volume in phase space. In the case of the 
microcanonical equilibrium this condition was met by the state that is maximally 
envariant, namely envariant under all unitary maps. Now, following Boltzmann’s line of 
thought we identify the canonical equilibrium by the configuration of the heat reservoir 
93 for which the maximal number of energy eigenvalues are occupied. Under the 
constraints, 

m m 

7 ij = N and Es — tk = A? e f (19) 

j = 1 i =1 

this problem can be solved by variational calculus. One obtains 

nj — n exp (Xef) , (20) 

which is the celebrated Boltzmann-Gibbs formula. Notice that Eq. (20) is the number 
of states in the heat reservoir 93 with energy ef for © and IB being in thermodynamic, 
canonical equilibrium. In this treatment temperature merely enters through the 
Lagrangian multiplier A. 

What remains to be shown is that A, indeed, characterizes the unique temperature of 
the system of interest, ©. To this end, imagine that the total system S can be separated 
into two small systems ©i and ©2 of comparable size, and the thermal reservoir, 93. It 
is then easy to see that the total number of accessible states 91 (e*,) does not significantly 
change in comparison to the previous case. In particular, in the limit of an infinitely 
large heat bath 93 the total number of accessible states for IB is still given by Eq. (18). 
In addition, it can be shown that the resulting value of the Lagarange multiplier, A, is 
unique [17]. Hence, we can formulate a statement of the zeroth law of thermodynamics 
from envariance - namely, two systems ©1 and © 2 , that are in equilibrium with a 
large heat bath 93, are also in equilibrium with each other, and they have the same 
temperature corresponding to the unique value of A. 

The present discussion is exact, up to the approximation with the Stirling’s formula, 
and only relies on the fact that the total system S is in a microcanonical equilibrium 
as defined in terms of envariance (6). The final derivation of the Boltzmann-Gibbs 
formula (20), however, requires additional thermodynamic conditions. In the case of the 
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microcanonical equilibrium we replaced conventional arguments by maximal envariance, 
whereas for the canonical state we required the maximal number of energy levels of the 
heat reservoir to be “occupied”. 

5. Concluding remarks 

In thermodynamics equilibrium states are characterized by sets of canonical variables 
[12]. For instance, for a closed system of given volume and given number of particles the 
internal energy, E, the volume, V, and the number of particles, N, are sufficient to fully 
describe all thermodynamic properties. These canonical variables are “macroscopic” 
in the sense that they are commonly interpreted as averages over many “permissible 
microstates”. From a “macroscopic” point of view it appears plausible that there may 
exist a vast number of microstates consistent with specific values of E, V, and N. In 
classical statistical mechanics [12] one then assumes ergodicity, i.e., that the individual 
system undergoes rapid, random transitions among these microstates, which is an excuse 
to ignore microstates and deal solely with the macrostate. This assumption is commonly 
formulated as the central postulate of statistical mechanics stating that [12] 

a macroscopic system samples every permissible quantum state with equal 

probability. 

Although this postulate appears to be plausible, and has proven to be powerful, it also 
poses several unanswered questions. For instance, it neither specifies the dynamics of 
the transitions nor does it specify the notion of probability that should be applied. 

Nevertheless, such formulations of the fundamental postulate have been used also in 
quantum analyses, e.g., to prove “Canonical Typicality” [15], or to study the interplay 
of “entanglement and the foundations of statistical mechanics” [14,18]. However, all 
these previous treatments had to rely to a certain extent on classical strategies that 
involved ensembles and in that context employed notions of probability and randomness. 
For instance, the starting point of the analysis in Ref. [15] is a set pure states, which 
is equipped with the Haar measure, i.e., an ad hoc selected probability distribution. 
Generally, however, probabilistic theories in physics are notoriously hard to formulate 
with mathematical rigor [11], Even two centuries after Laplace’s first description of the 
term “probability” [19] a conceptually universal approach to physical randomness - e.g., 
the relation between Kolmogorov-like measure theory and the underlying physical states 
- appears to be lacking within the classical context. Instead, several different notions 
of probability, such as subjective, objective, dispositional, frequentist etc. [20,21], have 
been developed, all based on non-equivalent mathematical concepts. 

The major achievement of the present analysis is a conceptually novel approach 
to the foundations to statistical mechanics. By exclusively relying on envariance - 
entanglement assisted invariance - we have identified quantum states that represent 
microcanonical equilibrium. The canonical equilibrium can then be described by simple 
counting arguments of degenerate energy states in a heat bath. Such counting is 
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reminiscent of classical statistical physics, but now we are counting not members of an 
ensemble, but provably equiprobable states. Our treatment is free of any mathematical 
or conceptual caveats. At no point ambiguous notions such as probability, ensemble, 
randomness, indifference, etc. have been employed. 
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Appendix A. Alternative treatment utilizing entropy 


In the preceding analysis we re-formulated the foundations of statistical mechanics in 
terms of envariance. This appendix is dedicated to an alternative treatment of (12), 
which involves the definition of the Boltzmann entropy FL. It is defined as the logarithm 
of the number of states, 

H = k B In (91) (A.l) 

where the parameter k b denotes the Boltzmann constant according to convention [12]. 
With this definition Eq. (12) can be re-written to read, 


«(«>) = exp (in 

= exp (k B ] FL^{E s - c fc ) - k^HsiEs)) . 


(A.2) 


We continue our discussion by inspecting each term in the exponent of Eq. (A.2) 
separately. The entropy of the total system S can be separated into a sum of the 
entropies of the subsystem of interest, (5, and the heat reservoir, 03, since we assume © 
and 03 to be uncorrelated. Accordingly, we have 


n s (E s ) = k B ln(m s (E s )) 

= k B ]n(Vle(E 6 )-Vl 8 (E tB )) (A.3) 

= U 6 (E 6 ) + / H<&{Et &), 

where Eq and E <g denote the internal energies of subsystem and heat reservoir, 
respectively. 

Now, we further assume that the internal energy of the heat reservoir is much 
larger than the one of the subsystem of interest, £# Eq — e fc , which is justified for 
dim(@) <C dim(OS) (see also Eq. (16)). Then, the first term in the exponent of Eq. (A.2) 
can be written as 
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where we identified E s — E & = E<$ and expanded the entropy of the heat reservoir 
around its average value. 

Substituting Eqs. (A.3) and (A.4) into Eq. (A.2) and making the usual 
thermodynamic identifications, &H^/dE^ = 1/T and Eg = Eq —T'Hq(Eq), we obtain 
the canonical state, 

01 (e fc ) = exp (^-^j . (A.5) 

As before, the crucial steps in the derivation are identifying Eq. (6) as microcanonical 
state and considering the fraction of states in Eq. (12). In this treatment, however, 
temperature is found to be a characteristic property of the heat bath. 

We emphasize again that both treatments are pure expressions of envariance and 
are free of ambiguous notions such as probability, randomness, indifference, etc. The 
main difference between the two presented derivations of the canonical equilibrium are: 
in the first treatment we obtained the canonical state for the heat reservoir 23, but we 
had to employ the additional thermodynamic assumption of maximally many energy 
levels being occupied; the second approach yielded the canonical state for © with no 
further information about 23, and we had to introduce the concept of entropy (A.l). 


Canonical equilibrium with system-bath correlations 


Finally, we briefly analyze the situation of non-vanishing coupling Hamiltonians, fig,®- 
For the sake of simplicity we still assume that we can neglect its contribution to the 
internal energies, but we relax the assumption of |sfc) being a product state (10). 
Generally, we have 

I s k ) = „(fc) \s u (k)) ® | b v (k)) , (A.6) 

V 

where | s u (k)) and \b u (k)) are again energy eigenstates in (5 and 23, respectively. This 
slight generalization allows for correlated and entangled 6 and 23. 

As before we are interested in the fraction of states available to S under the 
condition that the total internal energy, Eg, is given and constant (12). In contrast 
to the previous case (A.3), however, the total entropy FLs(Es) can no longer be written 
as a simple sum, since subsystem and heat reservoir are correlated. Instead we have 

H S (E S ) = He(Ee) + Hv(E v ) + k B In , (A.7) 


Note that the Shannon information of a maximally mixed states (6) is identical to the 
logarithm of the dimension of the Hilbert space, i.e., here the number of microstates [22], 
Thus, we can write with the reduced density matrices for the total system p$, the 
subsystem of interest pg, and the heat reservoir p®, 


D(p s \\pe ® p») = In 


/ 91s(Fg + E<s) \ 

v m^Eeym^j) ’ 


(A.8) 


where D(p s \\p 6 <g) p®) denotes the quantum mutual information [23]. 
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Accordingly, the generalized canonical state accounting for correlation between 
subsystem and heat reservoir becomes, 

9T(|s fc » = exp (-- D(p s \\pe ® P<b)^J • (A.9) 

Equation (A.9) generalizes the standard canonical state (A. 5) in a natural and intuitive 
way as correlations between subsystem of interest and heat bath are quantified by the 
quantum mutual information (A. 8). 
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